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Abstract 

Within our fission model, the Coulomb and proximity potential model (CPPM) cluster 

formation probabilities are calculated for different clusters ranging from carbon to silicon for the 

parents in the trans-tin and trans- lead regions. It is found that in trans-tin region the C, ^ O, 

TS[e and Mg clusters have maximum cluster formation probability and lowest half lives as 

compared to other clusters. In trans-lead region the ^^C, ^^'^°0, ^^F, ^^■^'^e, ^^'^°Mg and ^^Si 

clusters have the maximum cluster formation probability and minimum half life, which show that 

alpha like clusters are most probable for emission from trans-tin region while non-alpha clusters 

are probable from trans-lead region. These results stress the role of neutron proton symmetry and 

asymmetry of daughter nuclei in these two cases. 
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1, Introduction 

Cluster radioactivity which is the spontaneous emission of particles heavier than alpha 
particle can be treated as an intermediate process between alpha decay and spontaneous fission. 
This phenomenon was first predicted by Sandulescu et al. [1] in 1980 on the basis of quantum 
mechanical fragmentation theory (QMFT) [2]. After four years Rose and Jones [3] confirmed 
this phenomenon by the emission of C emission from Ra. Generally there are two 



approaches in explaining the cluster decay process: (1) cluster model and (2) fission model. In 
cluster model [4-6] the cluster is assumed to be preformed in the parent nuclei before it 
penetrates the nuclear interacting barrier. In fission model [7-9] the nucleus deforms 
continuously as it penetrates the nuclear barrier and reaches the scission configuration after 
running down the Coulomb barrier. In the cluster decay studies cluster-like shapes are preferred 
for very high asymmetric mass splitting and fissioning shapes are most suitable for less 
asymmetric and symmetric mass splitting. 

Clustering is a very general phenomenon, which appears in atomic, nuclear, sub nuclear, 
and the cosmic worlds [10]. A direct indication of clustering is the emission of alpha particle or 
heavier clusters from heavy nuclei. The properties of light nuclei with cluster structures have 
been well studied with cluster model. Clustering properties of nuclei would be helpfiil to obtain a 
systematic understanding of both the stable and exotic nuclei. Examples of such successful 
frameworks are the methods of fermionic molecular dynamics [11] and antisymmetrised 
molecular dynamics [12] which describe the structural properties of several nuclei, and their 
excited states, in the lighter mass region [13-16]. 

In nuclear dynamics, as seen in light stable nuclei, clustering is one of the essential 
features and various cluster structures have been known even in the low-energy region. Also, in 
the physics of unstable nuclei, clustering features comprise one of the central subjects. It is 
already weU known [12, 13, 17-20] that clustering structures appear in the ground states of 
ordinary light nuclei with A^= Z or in their neighborhood. Even though this phenomenon has now 
been studied for a long time, many things remain to be learned and recently some experimental 
activities focused on this subject [21] were reported. 



The present work aims to study the possibility of clustering in heavy nuclei by computing 
the cluster formation probability for the entire experimentally observed cluster decays from 
carbon to silicon taking the Coulomb and proximity potential as interacting barrier for the post- 
session region and simple power law interpolation for overlap region. We have also computed 
the cluster formation probability for different C, O, Ne and Mg clusters from '^ Ba, ^^ ' ^ Ce, 
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' Nd and ' Sm parents in the trans-tin region. The preformation probability can be 
calculated within a fission model as a penetrability of the internal part of the barrier, which 
corresponds to still overlapping fragments [22-24]. This part, however, is difficult to deal 
microscopically from current nuclear models due to the complexity of the nuclear many-body 
problem. 

In our previous works, Santhosh et al. [25-30] proposed the Coulomb and proximity 
potential model to study cluster radioactivity in various proton-rich parents with Z = 56-64 and 
N = 56-72, leading to doubly magic '°°Sn. We have studied the cold valleys in the radioactive 
decay of ^"^^"^^"^Cf isotopes [31] and the computed alpha decay half-lives closely agree with the 
experimental data. A semi-empirical model [32] is introduced for describing cluster radioactivity 
and this model has been applied to study the alpha decay of parents with Z = 85-102. The 
competition between the alpha decay and spontaneous fission half-lives of even-even nuclei with 
Z = 100-122 is studied [33] using the present model and the phenomenological formula for 
spontaneous fission. Recently we have proposed a semi empirical formula [34] for determining 
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the spontaneous fission half lives which is applied to the mass region from Th to 114. A 
brief description of cluster formation probability is given in section 2, the details of the Coulomb 
and proximity potential model are given in section 3 and the results and discussions are given in 
section 4. 



2, Cluster formation Probability 

The first model used to compute preformation probability was proposed by Blendowske, 
Flessbach and Walliser [4, 35, 36] on the basis of microscopic wave functions. For the case of 
C emission from Ra, the authors proposed an empirical formula for computing the cluster 
formation probability and is given as, 

S = {SJ'^-'''' (1) 

here S^ is the a-particle preformation probability and the values of S are different for even and 
odd nuclei. For even-even parent 5^''^" = 0.0063 and for odd-A parent nuclei, S"''' = 0.0032 [35]. 
Recently Dong et al. [37] computed the preformation probability using the eqn. 1 with the value 
of 5„ = 0.0338 for even-even parent and 5„= 0.0262 for odd-A parents. We have computed the 
cluster formation probability for C cluster emission from even-even parents by using the same 
relation and is found that ^'gig^jjo^j,,^ = 2.91x10 ,5^^^^=4.22x10 . 

Gupta and co workers [5] considered that in nuclear decay the formation of two 
fragments (the cluster and daughter nucleus) in their ground state with probability Po, defined the 
as a quantum mechanical probability of finding the fragments Ai and A2 (with fixed charges Zi 
and Z2, respectively) at a point R of the relative motion, the cluster formation probability is given 
as 

n(4)=K(4)fV^M^^ (2) 

where ^t^r; represents the mass parameter. The quantum number ft) = 0,1,2, counts the 

vibrational states y/^Z^ in the potential V. 



Within fission model Poenaru et al [38] defined the preformation probability, S as the 
penetrability through the internal part (overlap region) of the barrier: 

S = oM-KJ (3) 

Where 

K.=-\ ^2B{R){E{R)-Q)dR (4) 

Li which B{R) is the nuclear inertia, E{R) is the deformation energy and Ra the internal turning 
point, defined by E{R^ ) = , the energy released. 

There are many empirical relations put forward for computing the cluster formation 
probability, P^ based on different models. These are 

- log 10 Pq = 0.363 ^c + 3.044 for A^ < 30 Malik - Gupta [5] 

- logio Pq = 0.309 Ac + 2.214 Sandulescu et al [39] (5) 

- log 10 Pq = 0.733 Ac - 0.733 for even A Blendowske - Walliser [4] 

For Ac < 30 the three relations display the straight line description. For ^^ > 30 , only the 

calculations of Malik and Gupta are available for heavier cluster emissions and it is found that 
the cluster formation probability increases, rather than decreases, and becomes almost constant 
for ^^ > 40 . We would like to point out that in our earlier works [40] we have computed the 
cluster formation probability for all the clusters ranging from He to Ne for even-even parents 
with Z = 100-112 and found that P^ decreases with increase in mass number of cluster up to 

Ac =20 and then remains a constant. It is to be noted that different models [4, 5, 37-39] give 
different values for cluster formation probability but the computed T1/2 values agrees with 
experimental data which imply that Pg depends on models. 



3, The Coulomb and Proximity Potential Model 

The interacting potential barrier for a parent nucleus exhibiting cluster decay is given by 
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(6) 



Here Zi and Z2 are the atomic numbers of daughter and emitted cluster, 'r' is the distance 
between fragment centers, 'z' is the distance between the near surfaces of the fragments, i. the 
angular momentum, ju the reduced mass and Vp is the proximity potential is given by Blocki et 
al. [41] 



Vp{z) = A7iYb 
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With the nuclear surface tension coefficient. 



y = 0.9511 {1-1.1^26 {N-Zy I A'^ MeV/fm 



(8) 



Where N, Z and A represent neutron, proton and mass number of parent, O , represent the 
universal the proximity potential is given as [42] 



^{s) = -4A\e-'"''' 



for s>l. 9475 



(9) 



0(f) = -1 .7817 + 0.9270 s + 0.0169 s^ - 0.05148 s^ for < s < 1 .9475 



(10) 



With s = z/b, where the width (diffuseness) of the nuclear surface b ~1 and Siissmann central 
radii Ci of fragments related to sharp radii Ri is 
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For Ri we use semi empirical formula in terms of mass number Ai as [41] 



i?, =1.28X -0.76 + 0. 8X 
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(12) 



The potential for the internal part (overlap region) of the barrier is given as 



V = a,{L-L,y 



for z < 



(13) 



here L = z + 2C[ + 2C2 and Lq = 2C , the diameter of the parent nuclei. The constants «„ and n 

are determined by the smooth matching of the two potentials at the touching point. 
Using one dimensional WKB approximation, the barrier penetrability P is given as 



P = exp -|jV2//(^-0Jz 



(14) 



Here the mass parameter is replaced by reduced mass ju = mA^A2 / A , where m is the nucleon 
mass and Ai, A2 are the mass numbers of daughter and emitted cluster respectively. The turning 
points "a" and "b" are determined from the equation V(a) = V(b) = Q . The above integral can be 
evaluated numerically or analytically, and the half life time is given by 
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Where, u = 



^ co\ f2E.^ 



\2n J 



\ h J 



represent the number of assaults on the barrier per second and 'k the 



decay constant. Ey, the empirical zero point vibration energy is given as [8] 



E^ -e^ 0.056 + 0.039exp 
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for A2 > 4 



(16) 



For alpha decay, A2 = 4 and the empirical zero point vibration energy becomes, E^ = 0.095 Q 



The cluster formation probability, S can be calculated within the fission model as the 
penetrability of the internal part (overlap region) of the barrier given as 

S = exp(-K) (17) 

Where 

K^^]^2ju(V-Q)dz (18) 

a 

Here, a is the inner turning point and is defined as V(a) = Q and z = represents the touching 
configuration. 

4. Results and Discussion 

We have computed the cluster formation probability for the clusters ranging from carbon 

to silicon from the parent nuclei in trans-tin and trans-lead region using eqn. 18. Fig. 1 represents 

the plot for computed logio (S) versus A2, mass number of the fragment, for the emission of 

different carbon clusters from ^ ' ^ Ba, oxygen clusters from ^ '^ Ce, neon clusters from 

' Nd and magnesium clusters from ' Sm parents respectively. The Q values are 

computed using the experimental binding energies of Audi et al [43] and some masses are taken 

from the table of KTUY [44]. It is clear from the plot that the maximum value for cluster 

formation probability is obtained for ^ C, O, T^e and Mg clusters from the corresponding 

parents. One of the interesting facts is that the computed Q-values are high for these decays as 

compared to other reactions. As considered the neutron proton asymmetry distribution in light 

N 
clusters by the — ratio of the matter at the centre of the nuclei, for N = Z matter means that the 
Z 

neutron and proton distribution is symmetrical and only the alpha particle like matter exist at the 

centre of the nuclei. These nuclei are found to be in the neutron proton ratio equal to 1 and 



having greater cluster formation probability from the parent nucleus. Therefore we would like to 
point out that the clusters C, O, Ts[e and Mg form alpha like cluster configurations in the 
parent nucleus. 

Arumugam et al. [45] studied the possibility of clustering in light clusters with the 
increase or decrease of neutron numbers using relativistic mean field theory and examined the 
internal structure (density distribution) of various Be, C, O, Ne, Mg and Si isotopes. The authors 
found that the density distribution of Be indicates a a + 2p structure since there is more 

Q 

configuration of nuclear matter in the central region, the Be nuclei can easily be thought of as an 

o 

a-a cluster configuration. As the neutron concentration increases for Be nuclei, heavier than Be, 
the clustering in proton matter remains almost undisturbed, i.e., the core having alpha-alpha 

o 

cluster persists with the addition of neutron to the alpha-alpha nucleus Be, i.e. these isotopes 
have an a + a + xn structure, with alpha-alpha as the core. The remaining x neutrons are rather 
sparsely distributed around the core and constitutes a halo for these heavier isotopes of Be. For 
the case of carbon isotopes Freer et al. [18] and Horiuchi et al. [19] predicted the 3a-chain 
structure of '^C isotope. In the case of oxygen isotopes the ^''O is a doubly magic nucleus and it 
is spherical in shape. Zhang et al. [46] studied the structure of ^''O isotope and found that it forms 
a square, rhombic or a kite structure with alpha nucleus at its comer. The RMF calculations 
predict the ground state of O as having four alpha particles in the form of a regular tetrahedron 
in three-dimensional configuration and in two-dimensional configuration they form a kite 
structure, (see Fig. 6(b) of Ref. [45]) at deformation P = 0.95 and it yields the linear 4a-chain 
structure at highest deformation (P = 3.79). In 1994 Zhang et al [47] predicted ground state of 
TSle to be 5a-trigonal bipyramid and later the same structure was confirmed in many other 
calculations [4, 11]. Most of the calculations [47, 48] predict the ground state of Mg as a 



«+"'0 + «: structure with P~0.6 and having triaxial shape [49]. hi relativistic mean field 
mechanism the Mg isotope has the C+ C structure, though with difi'erent orientations - 
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prolate (central bispheroid) and oblate (trigonal biprism) and each C isotope is the combination 

98 

of three alpha particles. Many authors [16, 20, 47] predicts that the Si isotopes have the 
'^C + a+'^C trigonal biprism shape. From these observations we can see that the heavy clusters 
(from Be to Si) are formed as the combination of alpha particles and additional neutrons ( 
na+xn). 

Fig. 2 represents the plot for logio (S) versus A2, mass number of various carbon clusters 
from Ra, Fr, Ac and Th parents. It is evident from the plot that a peak in cluster formation 
probability is obtained at ^"^C nucleus. We would like to point out that the '"^C cluster emission 
from different Ra, Fr, Ac and Th parent nuclei has already been detected experimentally [50-55]. 
So the '"^C cluster preforms in these parents as compared to other carbon clusters, we would like 
to point out that the C cluster form (3o; + Iri) structure. 

Li Fig. 2 we have found a second peak at ^^C cluster but the ^''C cluster emissions are not 
observed experimentally. We have computed the half lives for ' ^^C cluster emissions from 
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Ra and Ac parents using Coulomb and proximity potential model. It is found that the 
lowest half lives is obtained for the ^ C cluster emissions from these parents. {For e.g. in the case 
of ^^^Ra, logio(7;/2) = 16.74 (for '^C) and log 10 (T;/^) = 27.45 (for ^^C); in ^^"^Ra, 

logio(7;/2) = 22.55 (for ^^C) and logio(7;/2) = 34.44 (for ^"C); in ^^^Ac, logio(7;/2) = 18.03 (for 

'"^C) and log 10 (ri/2) = 30.92 (for ^^C)} i.e. the '"^C cluster is most probable for emission, which 

preforms faster in these parents as compared to the C cluster. We would like to point out that 
many authors [56-58] predicts '"^C cluster emissions from these parents and our values are in 
agreement with their predictions. 
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The preformation probability of different oxygen isotopes from ' Th and " F 
clusters from Pa are shown in figure 3. It is obvious from the plot that in the case of Th 

90 

parent the maximum value of cluster formation probability is obtained for O and in the case of 
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Pa isotope the maximum value is obtained for F cluster, i.e., the cluster O forms (5a + 4n) 
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structures in the Th parent nucleus and F forms (4a + lp + 6n) structures in Pa isotope. It 
is obvious from the Figs. 2-3 that the cluster formation probability is greatest for C, O and 
O clusters from Th parent. The experimentally observed cluster decay for C and ^ O 
clusters has the same half life (experimental lower limit for log^ Tf,\^'{s) = 15.3 [55]). As we 
compare the preformation probability for these isotopes it is found that the cluster formation 
probability is maximum for '^C cluster (for '^C, 5 = 1.98x10^; for ^^0,5 = 2.27x10 'and for 
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0,5 = 2.98x10 ) . We have also computed the half lives for these cluster emissions using 
Coulomb and proximity potential model [32] and is found that the lowest half life time is for ^ C 
cluster emission (logjg 7^/2 (■5') = 18.74) as compared to '^0 (logjg T^,2(s) = 19.29) and ^"^O cluster 

emissions (log JO Ty2(s) = 25.26) . i.e., the C cluster preforms faster in Th parent as compared 

to other clusters. 

Fig. 4 represents the plot for logio(S) versus A2 for the emission of different neon clusters 
from various U, Th and Pa parents. It is clear from the plot that two peaks are present at A2 = 24 
and 26. In the case of Th and ' U parents the ' TMe clusters are experimentally detected 
[59-62] but TMe has the maximum preformation probability compared to Ne cluster. Again the 
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parents Th, Pa, ' ' U have the maximum cluster formation probability for Ne cluster 
as compared to T^e cluster, i.e these clusters are formed very fast in these parents as compared 
to other clusters. The Ne cluster has the (5a + 4n) structure and iS[e cluster has the (5a + 6n) 
structure. 
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For the case of U parent, the Ne and Ne cluster decays are experimentally observed 
with same half life (experimental lower limit for logjo Ty2(s) = 18.2 [55]). It is evident from the 

Fig. 4 that the Ne cluster has the greater cluster formation probability. We have computed the 
half lives for these decays using Coulomb and proximity potential model [32]; we found that the 
lower half life time is for Ne cluster emission (logjg Ty2(s) = 22.37) as compared to the Ne 
cluster emission (log 10 Ty2(s) = 22.60) . We would like to point out that the ^"^Ne cluster is formed 
faster than the other clusters from these parents and it has a (5a + 4n) combination. 

Fig. 5 represents the plot for logio(S) versus A2 for the emission of ^''"^^Mg clusters from 
various " U, Np, and ' Pu parents. It is clear form the plot that the highest peaks are 
found at ' Mg, the experimentally observed cluster radioactivity. In the case of " U, Pu 
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the maximum preformation probability is obtained for Mg cluster. For U parent ' Mg 
clusters have same half life (log 10 7;^f '(5) = 27.58 [63]) and in the case of ^^^Pu parent ^^'^°Mg 
clusters have same half life (logjo rif2'''(5) = 25.70 [64]) but the highest peaks in cluster 
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formation probability is obtained for Mg as compared to Mg cluster, i.e. the Mg cluster has 

9^7 ^0 

the greater chances for emission from these parents. In Np isotope the Mg cluster has the 
highest cluster formation probability. The Mg cluster forms (6a + 4n) combination while Mg 
cluster forms (6a + 6n) structure. 
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In the case of ' U parents the Ne and Mg cluster decays are experimentally 
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observed and it is clear that the lowest half life is obtained for Ne cluster (for U, 
log 10 T^iff (5) = 21 .08 [65] and for ^"u, logio Tiff (5) = 24.83 [66]). We would like to point out 

that our computation [32] also gives the same results (for UJogm Tyi2(s) = 20.72 and for U, 

log 10 Ty2(s) = 24.15 ), i.e. the Ne cluster is preformed faster in these parents as compared to 



other clusters. For the case of U parent the cluster radioactivity is observed for ' TSIe cluster 
with same half life (logioTiff (5) = 25.9)[60, 61] and for ^Vg cluster with 
logio 7lf2'"(5) = 27.54 [62], but the cluster formation probability is maximum for TS[e cluster 
compared to other decays. We have computed the half life for these cluster emissions and it is 
found that the lowest half life time is for ^Ne cluster emission (log jq Tj^jC^) = 25.88) as 

compared to ^"^Ne cluster decay (log jq Tj^j (5) = 27.39). hi the case of ^^''U isotope the ^"^e 
(logio^iff (5) = 25.9)[62] and ^^'^°Mg (log 10 Tiff (5) = 27.58) [63] cluster decays are 

experimentally observed, i.e the Ts[e cluster emission has the lower half life. As we compared 
the cluster formation probability of these clusters it is found that the Ts[e cluster has the larger 
preformation probability value. So we would like to point out that this cluster preforms first in 
the parent nucleus. 

The plot for logio(S) versus A2 for the emission of ^^"^''Si clusters from ^"^^Pu, ^"^^Am and 
Cm parents are shown in Fig. 6. It is obvious from the plot that the Si cluster has the 
maximum cluster formation probability in these parents and it has a {la + 6n) structure. 

hi the present study the maximum cluster formation probability is obtained for those 
decays which are experimentally observed. In brief our study reveals that ^"^C, ^^'^'^O, ^^F, ^'^'^'TSIe, 
' Mg and Si clusters have the maximum cluster formation probability in trans- lead region 
and ^ C, O, Ts[e and Mg clusters have the maximum preformation probability in trans-tin 
region and are due to the double shell closure of ^°^Pb and '°°Sn daughter respectively. Alpha 
like structures are probable for emission (maximum cluster formation probability) from trans-tin 
region and non-alpha like structures are probable for emission from trans-lead region which 
shows the role of neutron proton symmetry and asymmetry of daughter nuclei in these two cases. 
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Fig, 1, Plot for logio (S) versus A2, the mass number of various C, O, Ne and Mg clusters from 
Ba, Ce, Nd and Sm parents. 
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Fig, 2, Plot for logio (S) versus A2, the mass number of various carbon clusters from Ra, Fr, Ac 
and Th parents. 
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Fig, 3, Plot for logio (S) versus A2, the mass number of various oxygen clusters from ' Th 

9 1 9 S 9^1 

and " F clusters from Pa parents. 
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Fig, 4, Plot for logio (S) versus A2, the mass number of various neon clusters from U, Th and Pa 
parents. 
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Fig, 5. Plot for logio (S) versus A2, the mass number of various magnesium clusters from U, Np, 
and Pu parents. 
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Fig, 6, Plot for logio (S) versus A2, the mass number of various silicon clusters from ^"^"^Pu, ^"^^Am 
and Cm parents. 



